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Abstract 

A mechanism for relationship of solvent viscosity with reaction rate constant at 
enzyme action is suggested. It is based on fluctuations of electric field in enzyme 
active site produced by thermally equilibrium rocking (cranckshaft motion) of the 
rigid plane (in which the dipole moment ~ 3.6 D lies) of a favourably located and 
oriented peptide group (or may be a few of them). Thus the rocking of the plane 
leads to fluctuations of the electric field of the dipole moment. These fluctuations 
0^ \ can interact with the reaction coordinate because the latter in its turn has tran- 

sition dipole moment due to separation of charges at movement of the reacting 
£sj ■ system along it. The rocking of the plane of the peptide group is sensitive to the 

—^ , microviscosity of its environment in protein interior and the latter is a function 

of the solvent viscosity. Thus we obtain an additional factor of interrelationship 
for these characteristics with the reaction rate constant. We argue that due to the 
properties of the cranckshaft motion the frequency spectrum of the electric field 
fluctuations has a sharp resonance peak at some frequency and the corresponding 
Fourier mode can be approximated as oscillations. We employ a known result from 
the theory of thermally activated escape with periodic driving to obtain the reaction 
rate constant and argue that it yields reliable description of the preexponent where 
the dependence on solvent viscosity manifests itself. The suggested mechanism is 
shown to grasp the main feature of this dependence known from the experiment 
and satisfactorily yields the upper limit of the fractional index of a power in it. 



Key words: enzyme catalysis, Kramers' theory, thermally activated escape, 
periodic driving. 



Email address: sitnitsky@mail.knc.ru (A.E. Sitnitsky). 
Preprint submitted to 3 September 2008 



1 Introduction 



The functional dependence of the rate limiting stage k ca t for enzymatic and 
protein (ligand binding/rebinding) reactions on solvent viscosity r\ of the type 

k C at OC —; (1) 

rf 



where < (3 < 1 (usually (3 ~ 0.4 4 0.8) has been known for a long time [Tj, 
0, [3], 0], [5], [6], [7j, |B], [9], PS], P3]. More detailed studies revealed that 
in fact the fractional index of a power (3 is a function of cosolvent molecular 
weight M (i.e., the mass of a cosolvent molecule expessed in atomic units and 
measured in Daltons) [8] 

(3 = (3(M) (2) 



If one varies the solvent viscosity by large cosolvent molecules with high molec- 
ular weight that do not penetrate into enzyme then one obtains that the frac- 
tional exponent /3 — > 0, i.e., the reaction rate constant does not depend on 
solvent viscosity. With the decrease of cosolvent molecular weight the frac- 
tional exponent (3 increases. In the limit of hypothetical "ideal" cosolvent 
with infinitely small molecular weight (cosolvent molecules freely penetrate 
into enzyme and are distributed there homogeneously) it tends to a limit 
value Pmax = limM^o (3(M) ~ 0.79. The latter is neither experimental value 
nor a calculated one. It is an extrapolated number (see [8] for details). 



The functional dependence ([TJ also takes place for folding of proteins (see 
[T5] and refs therein). In our opinion enzyme catalysis and folding are quite 
different phenomena proceeding on different timescales (an enzyme turnover 
is typically 10 -4 4 10~ 3 s while "the time of folding varies from microsec- 
onds to hours" [E]). An enzymatic or protein reaction typically has a distinct 
rate limiting stage and can be perceived as an elementary step. Folding is a 
complicated process that involves huge number of elementary steps of com- 
mensurable importance. That is why we suppose that enzyme catalysis and 
folding involve different origin of the dependence ([TJ. In the present paper 
we deal only with enzymatic and protein reactions (i.e., those of bond break- 
ing/bond making in protein interior) and do not touch upon folding. For an 
unprejudiced observer folding seems to be an overworked issue in the litera- 
ture while physical aspects of enzyme action still remain in a deep shadow 
of their chemical counterparts for this phenomenon. Overwhwelming major- 
ity of researches from physical community perceive enzyme catalysis as some 
"chemistry" or "biology". That is why the aim to attract their attention to 
it as to a physical problem initiated by the the collection [JJ] and continued 



by the review article [UJ seems to remain as urgent as it was in the previous 
century. 

The attempts to explain the functional dependence ([TJ can be roughly di- 
vided into " phenomenological" and "theoretical". The former suggest that 
the fractional exponent /3 is the degree with which solvent viscosity is coupled 
with (frequency dependent friction) [15] or penetrates into (position depen- 
dent friction) [16], [T7] the protein interior. The latter try to derive it from the 
first principles [T5], [TS]. However Zwanzig model yields too small value for 
the fractional exponent (3 = 0.5 [T5]. Grote-Hynes theory [T5] gives that the 
rate dependence on solvent viscosity should be weaker than that predicted by 
Kramers' one (the latter yields k oc I/77 in the high friction limit [20] ). How- 
ever no explicit derivation of expression (pQ) from the Grote-Hynes theory has 
been achieved. As the authors of [S] conclude "there seems to be no general 
agreement yet about the origin of the fractional (3 value in Eq.l". The authors 
of [10] draw to a similar conclusion. In our opinion little has changed in this 
issue (as applied to enzyme catalysis only because there is certain progress 
in understanding of viscosity dependence for folding [13]) since the date of 
the cited papers. The aim of the present paper is to provide theoretical inter- 
pretation of the functional dependence ([TJ) and to "explain" the limit value 

Pmax ~ yj.lt). 

There seems to be a consensus among researhes in understanding that the 
dependence of an enzymatic reaction rate constant on solvent viscosity is me- 
diated by internal protein dynamics. This undersanding goes back to the so 
called transient strain model. The latter is based on the idea of overcoming 
the energy barrier of an enzymic reaction by structural fluctuations whose 
frequency is inversely proportional to the viscosity of the medium [21], [22], 
[23] . That is why any theory of the phenomenon should be a part of the main- 
stream of modern enzymology to study the role of dynamical contribution 
into enzyme catalysis (see the materials of a recent conference in the sub- 
ject issue of Phil. Trans. R. Soc. B (2006) 361). There are different sonorous 
names for such dynamical mechanism: "rate promoting vibration" (RPV) |24j . 
the "protein promoting modes" [25], [26], etc. In the present paper the name 
RPV is used as the most appropriate one for the concept under consideration 
that some conformational motion of vibrational character in protein is coupled 
somehow to the reaction coordinate. However it should be stressed that the 
author of the present paper input in this name absolutely different meaning 
than the authors of [21], [25], [2S] and other papers within the framework of 
this concept. We invoke to the idea that a dynamically unusual electric field 
in enzyme active site may play a key role for catalysis. This idea was put for- 
ward by Frohlich in his concept of coherent vibrations of protein giant dipole 
moment [2Z] , [2H] , Gavish and Werber in the hypothesis of charge fluctuations 
[TJ and Warshel in his concept of electrostatic fluctuations [22] , [30] (see also 
and refs. therein). In the present paper the name RPV means 



the following: a Fourier mode of the fluctuating electric field in the enzyme 
active site generated by protein dynamics [34] , [35] . 

Warshel and coauthors [31], [22], [33] argue the following statements. 1. A 
dynamical mechanism can contribute significantly into enzyme catalytic effi- 
ciency only if it leads to nonequilibrium (non-Boltzmann) distributions for the 
reaction coordinate produced by coherent oscillations in protein dynamics. In 
the opposite case of thermal equilibrium dynamical effects can lead to noth- 
ing more than some modest corrections in the preexponent. 2. Equilibrium 
protein dynamics can not lead to coherent oscillations coupled to the reaction 
coordinate. The item 2. from this list is doubtless. However the item 1. in 
our opinion is not so and an efficient dynamical mechanism can stem from 
thermally equilibrium fluctuations. Moreover even if the item 1. is true (i.e., 
a dynamical mechanism does not contribute into the catalytic efficiency) the 
corrections in the preexponent can be crucial for the dependence of the enzy- 
matic reaction rate constant on solvent viscosity because the latter manifests 
itself namely in the preexponent. We discuss a possibility that protein dynam- 
ics produces specific fluctuational influence on the reaction coordinate. This 
influence on the one hand is of thermally equilibrium origin and on the other 
hand is additional to those available for reactions in solution (i.e., the latter 
have no analogous counterpart in the thermal noise spectrum). Of course this 
fluctuational influence can not be coherent oscillations. However in our opinion 
coherent oscillations are not crucial to be the origin of the dynamical mecha- 
nism. Fluctuations of thermally equilibrium nature can play this role as well. 
Regretfully as will be argued below it seems rather difficult to treat such fluc- 
tuations in their natural form. That is why in the present paper we invoke to 
the fact that some Fourier mode from their frequency spectrum mimic coher- 
ent oscillations so much that can be considered in the first approximation as a 
steady harmonic vibration. We stress once more that it is merely a methodical 
trick to reduce the problem to elaborate theoretical technique rather than an 
indispensable assumption for the present approach. The question "how much 
can such vibration contribute into the reaction rate enhancement ?" is not the 
matter of the present paper and is touched upon rather briefly here. The aim 
of the paper is to show that this vibration (that is the RPV in our approach) 
enables us to interprete experimental data on solvent viscosity depenedence 
for enzymatic reactions. 

The paper is organized as follows. In Sec. 2 the relevant protein dynamics as 
the origin of the electric field fluctuations is introduced. In Sec. 3 the interac- 
tion of these fluctuations with the reaction coordinate is disscused. In Sec. 4 
the equations of motion are obtained and the reasons why the electric field 
fluctuations can be conceived as oscillations are argued. In Sec. 5 the influ- 
ence of these oscillations on the reaction rate is considered. In Sec. 6 numerical 
estimates are presented. In Sec. 7 the solvent viscosity dependence for the re- 
action is obtained. In Sec. 8 the results are discussed and the conclusions are 



summarised. In Appenix A some technical details are presented. In Appendix 
B some additional material is presented. 



2 Origin of the RPV 



A crucial question for the mainstream of modern enzymology to investigate 
the role of dynamical effects at enzyme catalysis is the following: how does 
the RPV that is typically on the picosecond time scale affect the catalytic act 
that is typically on the millisecond time scale (an enzyme turnover is usually 
~ 10~ 3 -7- 10 _4 s) (see [23] and refs. therein)? In our opinion the frequency 
of the RPV as itself is not an appropriate characteristic for this issue. As 
a matter of fact it is not of much significance whether the RPV is on the 
picosecond time scale or, e.g., on the nanosecond time scale. The relevant 
and the most important characteristic is the life time of vibrational motion 
in protein dynamics that produces the RPV. Namely the problem of survival 
of vibrational excitations on the time scale of enzyme turnover plagues many 
of the speculations about dynamical contribution into enzyme action and is a 
point of application for criticism by Warshel and coauthors [31] , [32] , [33] . It 
is argued below that this problem does not arise in the present approach. 

The most natural candidate for the source of electric field is a peptide group of 
protein backbone. The latter is known to have a rather large constant dipole 
moment p ~ 3.6 D that lies in its plane [36]. The dipole moment produces an 
electric field. Thermal fluctuations (rocking) of the rigid plane of the peptide 
group relative to its mean averaged position in the protein backbone lead to 
variation in time of this electric field, i.e., to the electric field fluctuations. As 
the latter are produced by thermally equilibrium fluctuations they exist on 
the whole duration time of the catalytic act. That is why there is no problem 
to match the electric field fluctuations with the process of catalysis. We are 
interested in the amplitude and spectral properties of these fluctuations in the 
enzyme active site at the place of the reaction coordinate. 

The rocking of the rigid plane of the peptide group (the so called "crankshaft- 
like" motion) due to degrees of freedom of torsional (dihedral) angles (pi and 
ipi-i (see Fig.l) was proposed on theoretical grounds from normal- mode anal- 
ysis [37J. It is supported by numerous NMR experiments and molecular dy- 
namics simulations of protein backbone [SB], [3H], [3D], [31], [32], [33]) [33], [32], 
p£6] , [3Z] , [3H] , [32] , [50] , [SI] . The crankshaft-like motion is comprehended now 
as a dominant type of motion for the protein backbone that "involves only 
a localized oscillation of the plane of the peptide group" [51] . The essence of 
this motion is the so called anticorrelated motion of the torsional angles (pi 



and ipi-i manifested itself in the requirement (see Fig.2) 

0/2 = Aip, = AV»i_i (3) 



In this case the plane of the peptide group rocks as a whole around some 
axis a that goes through the center of masses of the peptide group parallel to 
the bonds C % ~ x — C % and N l — C l a (see Fig.l). The moment of inertia of the 
peptide group relative to the axis a is known to be I ~ 7.34 • 10~ 39 g ■ cm 2 
[52] . However this value is given in [52] without a reference or calculation. 
That is why in Appendix A a brief estimate corroborating this value is given. 
Molecular dynamics simulations and NMR experimental data testify that the 
character of the correlation function for the crankshaft motion is decaying 
oscillations [ID], [51]. However they provide characteristics for them in a very 
wide range from the subpicosecond and picosecond time scale [51] to slower 
motions on a much larger time scale from tens of picoseconds to 100 ps and 
more [50]. This is presumably a particular manifistation for the cranckshaft 
motion of a general principle of hierarchical structure for a conformational 
potential in protein dynamics [3J, [5], [IT] when a group of local minima forms 
a smooth local minimum and so on. In fact knowing the actual values of the 
frequency for the oscillations and the characteristic time of their decay for the 
functionally important crankshaft motion is not indespesible for the purposes 
of the present paper. However in our opinion it is reasonable to assume that 
the frequency of oscillations of the plane of the peptide group as a whole for 
such motion should be at least order of magnitude less than those for high 
frequency in - plane motions such as, e.g., Amide I (~ 1600 cm^ 1 ). The choice 
of the frequency in the ujq ~ 100 cmT x m 10 13 s^ 1 range enables us to match 
it with the amplitudes of rocking of order of several degrees (see (1T21) below) 
in accordance with experimental data |47| . [ID] . 

It is natural to describe the cranckshaft motion by a Langevin equation. Such 
equations are frequently used in protein dynamics [53J, [53], [55]. The equation 
of the cranckshaft motion for the rigid plane of the peptide group in the above 
conditions is 

/^ + 7^ + /<*(«)=««) (4) 



where 7 is the friction coefficient and £(£) is the random torque with zero 
mean < £(£) = > and correlation function 

< £(0)£(*) >= 2k B Tj5(t) (5) 

ks is the Boltzman constant, T is the temperature. 



We consider the hydrodynamic friction. This "macroscopic" notion is known to 
work surprisingly well at the molecular level (see [56] for thorough discussion) . 
We model the peptide group by an oblate ellipsoid with halfaxes a, b and c 
(where a ~ c and a,c >> b that reflects the flat character of the peptide 
group) Fig. 3. For the rotation of the ellipsoid around x-axis (that is in our 
case actually the axis a for the peptide group introduced above) the friction 
coefficient is given by a formula [57] 



(7)x = 87r?7 abc 
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where r = 2(abc) 1 ^ 3 and i] is the viscosity of the ellipsoid environment. The 
required behavior is obtained if we have the condition of the underdamped 
motion 
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(7) 



Let us meke a numerical estimate of this requirement. Taking the linear size 
of the peptide group (distance for the atom O to the atom H) to be L p» 2.5 A 
we have c ~ L/2, a ~ c/2 and b ~ a/3. From (Q we obtain that at r\ « 1 cP 
(that is the viscosity of water at room temperature) 7^2- 10~ 27 CGS. Then 
at typical value ujq ~ 10 13 s _1 and the value I w 7.34 • lCT 39 g 



cm 



we have 

7/(iu>o) ~ 0.01, i.e., the requirement ([JJ is by far satisfied. Moreover it well 
holds even at the increase of 77 by an order of magnitude. 

At the requirement (J7_D we obtain (see e.g. [58j, [55] ) 



a(t) =< 0(O)0(t) >f 
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(9) 



Then the Fourier spectrum for the correlation function a(t) is 



a(cu) 
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(10) 



At our requirement ([7]) it has a sharp resonance peak at the frequency io>o- For 



the mean squared amplitude (msa) we have 

msa = — 5- (11; 



The latter means that the amplitude of the cranckshaft motion at room tem- 
perature and, e.g., typical value u ~ 10 13 s _1 is 



lk B T 



0.1^6° (12) 



that is Aifi = Aipi_i « 3°. This value sheds light on the origin of essentially 
vibrational character of the peptide group motion manifested itself in (JSJ) and 
(fTUl) . At such angles of rotational deviation of the peptide group from its mean 
averaged position the linear displacements of the atoms are ~ c-<f) max m 0.1 A. 
The latter value is much less than both the size of a solvent molecule and 
interatomic distances to neighbor fragments of protein structure. That is why 
the environment exerts rather weak friction for such type of the peptide group 
motion that is reflected in the requirement (JTj). Thus we conclude that the 
thermally equilibrium cranckshaft motion of the peptide group is of essentially 
vibrational character even in such condensed medium as protein interior. 



3 Essence of the RPV 



We consider the plane of the peptide group that undergoes thermally equi- 
librium rocking around its mean averaged position. The angle <j)(t) quantifies 
the random deviations of the plane. We choose the axis x in the direction of 
these deviations (see Fig. 4). Now we recall that there is the dipole moment of 
the peptide group with the absolute value p ~ 3.6 D laying in its plane. We 
choose the axis z in the direction of the dipole moment at the mean averaged 
position of the peptide group. Thus the dipole moment is a vector in our frame 

pit) = p sincj)(t) e x + p cos<f>(t) e z (13) 



Taking into account that | <f>(t) \« 1 we have 

p(t) w p 4>{t) e x + p e z (14) 



We assume a simplest geometry favorable for catalysis (generalizations are 
trivial but lead to more cumbersome formulas of not principle character) . In 



this geometry at a distance R from the dipole moment along the axis x 

R = Re x (15) 

there is a reaction coordinate q. The latter has a transition dipole moment d 
due to separation of charges with the value e at movement along it 

d = e(q-q )e x (16) 



Here go is the position of the minimum of the potential energy surface corre- 
sponding to initial reagent. Then the quasistationary fluctuating electric field 
of the peptide group dipole moment in the place of the reaction coordinate is 



3 (p(t) -R)R- R 2 p(t) 



m = ^— — ^g — — * ^s W)e. - e ~*l ( 17 ) 



where e is the dielectric constant. Denoting 
2ep 



9 



eR 3 



we obtain for the interaction of the peptide group dipole moment with the 
transition dipole moment of the reaction coordinate 

H int = -(d-E(t)) = -g(q - q )<j>(t) (19) 



This interaction reveals the essence of the RPV as fluctuating electric field 
affecting the reaction coordinate. 



4 Equations of motion 



We denote m the effective mass of the reaction coordinate q and V(q) the po- 
tential energy surface for it. We assume the most frequently used memory-free 
friction (Ohmic damping) for both the reaction coordinate and the peptide 
group cranckshaft motion. For the reaction coordinate it can be motivated 
as follows. Taking into account the memory friction (that leads to the gen- 
eralized Langevin equation, e.g., in the spirit of Grote-Hynes approach) is 
usually necessary for fast reactions (e.g., those with low activation barriers). 
Enzymatic reactions typically have a rather high barrier (even after its re- 
duction by catalytically active groups). Thus the limit of a very rapidly fluc- 
tuating force compared with the rate of the transition over the barrier (the 



reaction dynamics is very slow and the environment exerts its full frictional 
influence during barrier crossing) can be safely taken. In this case the gener- 
alized Langevin equation is well approximated by the ordinary one. For the 
high - frequency cranckshaft motion of the peptide group the approximation 
of Ohmic damping can not be logicaly motivated because the characteristics of 
the correlation function (inverse frequency of oscillations and decay time) are 
comparable with the characterisic time of surrounding rearrangements acting 
as thermal bath. Besides taking into account memory friction for peptide dy- 
namics can be implemented rather easily in the approximation of harmonic 
conformational potential [15] . [51], [5_H], [SO], [SI]- However it usually leads 
to complications of not principal character at room temperatures (well above 
the so called "glass" transition in protein dynamics) where enzymes normally 
work. In this temperature range conformational protein dynamics are satis- 
factorily described by ordinary Langevin equation. 

Then we have classical equations of motion are 

mq + V'(q) -g<p+vq = ((t) (20) 

70 + Iul<P - g(q - g ) + 7 = £(t) (21) 



where the coupling constant g is given by (ITS]) . In (1201) v is the friction co- 
efficient for translational motion along the reaction coordinate (modeled by 
effective Brounian particle with the caharacterisic linear size / in the media 
with viscosity rj and given by the Stokes formula v = Qirlr} (see next Sec.)) 
and ((t) is the random force characterized by 

< C(£) >= 0; < C(0)C(t) >= 2k B Tu5(t) (22) 



In (1211) 7 is the friction coefficient for rotational fluctuations of the peptide 
group given by (Q and £(£) is the random torque characterized by 

< m >= 0; < £(0)£(i) >= 2£; B T7<5(£) (23) 



Neglecting the "backward" influence of the reaction coordinate on the peptide 
group motion (term —g(q — qo) in (|2~T|) ) we obtain for the latter the previously 
considered equation of motion (j4j) and are left with a two noise problem 

mq + V'{q) + vq = C(t) + g<f>(t) (24) 



Here the internal (its intensity is related with friction coefficient v) white 
noise ((t) is characterized by (I2"2"j) and the external (it does not create friction 



10 



for the movement along the reaction coordinate q) oscillating noise <p(t) is 
characterized by 



< 0(t) >= 0; a(t) =< 0(O)0(£) >^ j^exp I -- ^ I cos (uj t) (25) 



k B T_ ( 7 | t 



The stochastic influence g<fi(t) has no counterpart for reactions in solution. 
It is unique for enzymatic reactions because it is produced by dynamics of 
protein structure (namely by a favourably located and oriented peptide group 
in the present model or may be by a few of them in more realistic cases). In 
solution solvent molecules also posess dipole moments and undergo thermal 
motion. However they do not form strictly determined structure enabling them 
to implement high-frequency motion of essentially vibrational character as in 
the case of peptide groups in proteins. 

The equations fl2^|) . ff22l . ([25]) belong to a class of the problems considered 
recently in [62] , [63] . Regretfully the formulas obtained there lead to very 
cumbersome manipulations in our case and this most natural way of formu- 
lating the problem has not led to representable results for oscillating noise 
(1231) yet. That is why we have to resort to simplifying assumptions. First of 
all we take into account that a reaction of bond breaking or bond making 
requires linear displacements of atoms of order of the bond length ~ 1 A 
that is comparable with the size of solvent molecules. Hence the solvent exerts 
rather strong friction for the movement of the system along the reaction coor- 
dinate. Thus we can restrict ourselves by the high friction limit, i.e., neglect 
the inertial term mq in (T2^|) . The requirement for the overdamped regime is 
vjivnio^ » 1 where the friction coefficient is given by the Stokes formula 
v = 6nlr] (see next Sec.) and the frequency uj^ characterizes the shape of the 



potential barrier at its top uj^ = J\ V"(qb) | /m. At I pa 1 A and r\ ~ 1 cP 
(that is the viscosity of water at room temperature) we have the estimate 
v pa 2 • 1CT 9 CGS. Then at typical values u b ~ 10 3 cmT x pa 10 14 s _1 and 
m pa 5 a.u.m. (a.u.m. pa 1.7 • 10~ 24 g) we obtain vjimio^) pa 2. Thus the 
requirement of overdamped regime does not hold well in pure solvent but be- 
comes satisfactory with the increase of viscosity. A generalization of the theory 
with taking into account the inertial term in the ordinary Langevin equation is 
indispensable for the description of reactions in a gas phase that proceed in the 
underdamped limit. Those in condensed media (in solution or in an enzyme) 
are generally believed to proceed typically in the overdamped regime and we 
also resort to this assumption for simplification of further analysis. We stress 
once more the reason why for the motion of the peptide group we apply the 
underdamped regime while for the movement along the reaction coordinate 
we invoke to the overdamped one. In the former case the linear displacements 
of the atoms are pa 0.1 A, i.e., much less than the size of the solvent molecules 
pa 1 A (see Sec. 2). In the latter case the linear displacements of atoms are 
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comparable with the size of the solvent molecules. That is why the friction is 
supposed to be sufficiently strong for the overdamped regime to be applicable. 

Most important of all we resort to a severe approximation based on the fol- 
lowing reasoning. Strictly speaking a Fourier mode of the random process 4>(t) 

$(u) = — f dt <f>(t)exp(-iwt) (26) 

2ir J 



is a random function of frequency with zero mean. However the mean squared 
amplitude of this mode is related to the correlation function a(t) via the 
Wiener-Khinchin theorem 

<| 0( w ) | 2 >= / dt a(t)exp(-iut) = 2na(uj) (27) 



where a(u) is given by (flOl) . We can introduce the effective mean amplitude 
at any frequency as the square root of (l2Tj) . In particular we can do it for the 
frequency u 



<| 0(w o ) \>eff= \J2-Ka(u ) (28) 

As was already mentioned after (fit)]) the spectrum of a{uo) at the requirement 
([7]) has a sharp resonant peak at the frequency ujq. It means that practically 
the whole power of the spectrum is in the Fourier mode <P{ujq) of the random 
process <p{t) at this frequency. That is why we can roughly approximate the 
random process (p(t) by a harmonic vibration with the amplitude given by 



<p(t) «<| 4>(u ) \> eff sin(u t) (29) 

From ([H]), (TTUT) with taking into account ([7]) we obtain 

0(«) 



\ 



V2k B T . 

stn{uJot) [SO) 



7^o 



This artificial replacement of the oscillating noise (that is a random process 
with the properties (I2"oD ) by oscillations (I5U1) simplifies the problem signifi- 
cantly and enables us to employ some known results from the Kramers' the- 
ory. As was stressed above it is a methodical trick enabling us to cast the 
problem into a tractable form rather than an indespensable assumption for 
our approach. 
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5 Effect of the RPV 

A standard tool to investigate dynamical effects in reaction rate is the Kramers' 
theory (for review see [20], [HI], [SS], [HE] and refs. therein). The latter is based 
on the ordinary Langevin equation of motion along the reaction coordinate 
q. In the Kramers' theory a chemical reaction is modeled as the escape of a 
Brownian particle with the mass m (the effective mass of the reaction coordi- 
nate, e.g., the reduced mass of a scissile bond) and linear size I in the media 
with viscosity rj (so that the friction coefficient is given by the Stokes formula 
v = 6ttIt]) and temperature T from the well of a metastable potential along the 
reaction coordinate. This potential is an electronically adiabatic ground state 
obtained in the Born-Oppenhaimer approximation by methods of quantum 
chemistry. It is considered as preliminary input information for the Kramers' 
theory and neither its origin nor its possible modification by catalysts is an 
issue of the present approach. 

For the potential V(q) we consider an arbitrary form with a barrier that has at 
the top of the latter the value V"(qb) ~ -mwj and at the bottom of the well 
the value V"(q a ) ~ rnu^. Also we include into consideration the oscillating 
force 

g<f>(t) = h sin(u>ot) (31) 

where 



h 



2ep 

7ip\ 



^" T (32) 



7^0 



As a result the ordinary Langevin equation of motion along the reaction co- 
ordinate q (1241) takes the form 

v-j- = -V'(q) + h sin(u t) + ((t) (33) 

LLC 



We introduce the dimensionless variables and parameters as follows 

Q -*■ q(Qb ~ q a ); t-^t J? ^o ->• 2 =fi ( 34 ) 



and denote 

h ksT 

mul(q b - q a ) ' mu^{q b - q a ) 2 
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where A is the dimensionless parameter characterizing the oscillating electric 
field strength and D is the dimensionless parameter characterizing the inten- 
sity of thermal fluctuations (i.e., the temperature of the heat bath) relative to 
the barrier height. Their ratio is the most important parameter of the model 



.4 
D 



H<ib - la) 2ep{qb - q a ) 



k B T 



eR 3 \ ju k B T 



V2 



(36) 



In dimentionless variables the Langevin equation is 
j t =-U\q) + Asin(Slt)+ X (t) 



(37) 



where U(q) is the dimensionless potential and the dimensionless noise x{t) — 

rwtd-^ c (£r) is characterized b y < *(*) >= ° and < x(o)x(*) >= 

2DS(t). The corresponding overdamped limit of the Fokker-Planck equation 
for the probability distribution function P(q,t) [SD] is 



dP(q,t) d 



dt 



— {[-[/'(g) + Asin(Qt)]P(q,t)} + D 



dq 



d 2 P(q,t) 
dq 2 



(3J 



In the absence of driving (A = 0) the escape rate is given by the famous 
Kramers' formula 



K 



LJaUb 
27T 



-exp 



(U(q b )-U(q a ))/D 



(39) 



where u a = Jll"(q a ) and u b = J\U"{q b )\. In the case of thermally activated 
escape with periodic driving (A ^ 0) one usually introduces the instantaneous 
rate constant T(t), the rate constant T averaged over the period of oscillations 
T = 27r/Q and is interested in the escape rate enhancement T/T K (see [H], 
[67] , [68] . [69] and refs. therein). Regretfully no workable formula for the gen- 
eral case of arbitrary modulation amplitude to noise intensity ratio A/D and 
arbitrary frequency Q is available at present. However the case of moderately 
weak to moderately strong modulation D « A « y/D is relevant for our 
problem. In this case a simple formula for the escape rate enhancement is 
known [68] 



-I 



2nAc(Q)(l + 0(A)) 



-.exp 



Ac (Q) 
D 



(1 + 0(A)) 



(40) 



where for the cubic (metastable) potential c (Q) = irfl/sh (jrfi) and for the 
quartic (bistable) potential c (fi) = JTc£l/sh(ir£l). This formula is sufficient 
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for the purposes of the present paper to understand the solvent viscosity de- 
pendence in the preexponent. However it should be stressed that it is invalid 
is one wish to evaluate reaction rate enhancement due to the suggested mech- 
anism. In this case the corrections 0(A) lead to deviation from the log-linear 
behavior predicted by fj40l) . Fortunatly the latter problem is not a matter of 
the present paper and (J4"0j) suffices for our needs. However we briefly return 
to the problem of evaluation of the reaction rate enhancement due to the 
suggested mechanism in the next Sec. and in Appendix B. 



6 Numerical estimates 



Before proceed further we should define the range of the parameters for our 
model. To be supported by evidence we take the numerals for a very typical 
and one of the most studied enzymatic reaction catalysed by Subtilisin. This 
enzyme belongs to serine proteases and brakes the bond between the atoms 
C and N in a peptide group of a substrate of protein nature. At physiological 
temperatures where enzymes normally work (ksT « 4.2 • 10~ 14 CGS) the 
corresponding noncatalysed reaction in solution typically has the rate constant 
k « 1 • 10~ 8 s~ x while the enzymatic reaction has the rate constant (that of 
the rate limiting step) k cat ps 5 • 10 1 s~ 4 [70]. Thus the total catalytic effect 
is approximately 10 12 . We have for the friction coefficient of the movement 
along the reaction coordinate v = 6ttIi] at I ~ 1 A and 77 « 1 cP (that is 
the viscosity of water at room temperature) the estimate v m 2 • 10~ 9 CGS. 
Then taking into account that we measure dimensionless time in the units 
of — ^ (see ( 1341) ) we have for the reaction rate constant of the noncatalysed 

b 

reaction in dimensionless form at m « 5 a.u.m. (a.u.m. w 1.7 ■ 10~ 24 g) and 
ujh ~ 10 14 s _1 the value V = -^W « 10~ 22 . We can evaluate from the Arrhenius 

factor of the Kramers' rate in the dimensionless form for the quartic (bistable) 
potential V(q) = -g 2 /2 + g 4 /4 [20] 



T K = (V2n)- 1 



■k) exp 



■1/(4D) 



(41) 



that in order to obtain such a reaction rate constant in the absence of the 
oscillating electric field (A = 0) we should take D w 422 1 fal0 w 5 ■ 1CT 3 . If 
we take the cubic (metastable) potential V(q) = q 2 /2 — q 3 /3 then we obtain 
D«3' 10~ 3 . Thus the typical noise intensity for enzymological problems is 

D ss 3 4- 5 • 1(T 3 (42) 



As earlier we assume that the dimensional frequency of the RPV is on the 
picosecond time scale and is approximately u rj 100 cmr 1 w 10 13 s^ 1 . Then 
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from (1341) we obtain at m 
frequency in our case is 

fi^-KT 1 « 1 



5 a.u.m. and uj b 



10 14 s x that the dimensionless 



(43) 



Thus we are actually in the low frequency regime or that of slow modulation 
Q << 1. In this case we have c(O) ~ 1 for both cubic (metastable) and 
quartic (bistable) potentials. That is why c (Q) will be droped further from 
the formula for the reaction rate enhancement. 



We take the value p = 3.6 D and typical values: for the partial charge 
e ps 0.3 a.u.c. (a.u.c. = 4.8-1CT 10 CGS), for the reaction potential q b — q a ps 1 A, 
for the dielectric constant in protein interior e ps 3.5, for the distance of the 
fluctuating peptide group from the reaction coordinate R ~ 3 A, for the fre- 
quency of the RPV ujq w 100 cm -1 ps 10 13 s _1 and for the friction coefficient 
of the peptide group (see Sec. 2) 7^2- 10~ 27 CGS. Then at room temper- 
ature we obtain from (|36|) that A/D ps 10. At such value of this parameter 
the formula (j4"0~]) predicts the reaction rate enhancement by 3 -j- 4 orders of 
magnitude (log A ps 3.4). This estimate requires comments. As was stressed 
above the formula (140]) overestimates the reaction rate enhancement. As a 
matter of fact the corrections O(A) lead to some more moderate growth than 
the log- linear one [57], [55] , [55] . However the results presented in these pa- 
pers testify that an appreciable deviation from the log-linear behavior begins 
at significantly higher values of modulation amplitude to noise intensity ratio 
A/D than A/D ps 10. That is why the reaction rate enhancement by ~ 3 -j- 4 
orders of magnitude at A/D ps 10 due to the suggested mechanism seems 
to be quite feasible. Here we restrain ourselves from further discussing this 
mostly important issue because the present paper is devoted to a particular 
effect manifesting itself in the preexponent rather than in the exponent of the 
formula (14*0]) . However in Appendix B we return to this problem and discuss 
it in more details. 



7 Solvent viscosity effect 



Returning to dimensional parameters, substituting (|36|) and (|39|) into (|4"0~]) and 
dropping c (O) ~ 1 (see previous Sec.) we obtain for the reaction rate constant 



r = r K A 



my/V"(q a )\V"(q b )\ 



1-KV 



exp 



V(q b ) ~ V(q a : 



eR 3 



l^ujQk B T 



\ A7rep(q b - q a )\ ^2 



-exp 



2ep(q b - q a ) 



eR 3 \ "yuj k B T 



V2 



x 



(44) 
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First we consider the hypothetical "ideal" cosolvent with infinitely small molec- 
ular weight (cosolvent molecules freely penetrate into enzyme and are dis- 
tributed there homogeneously). Taking into account that in this case we have 
both v oc r] [y = 6nln) and 7 oc 77 (see (JH])) we obtain 

r « ^75 ( 45 ) 

Thus the present approach yields the limit value j3 max = 0.75 for the "ideal" 
cosolvent that is rather close to the extrapolated limit value f3 max ~ 0.79 
from the paper [8]. The dependence on 77 in the second exponent is strongly 
overshadowed by the first exponent because typically we have 



V{%) ~ V(q a ) ^ ^ , 2ep(q b - q a ) 



V2 



kb-L en 6 \ bivluJokBi 



Then in the In of the reaction rate constant InT = const — 0.75 In n+B — b/^/rj 
the last term is negligible provided the requirement (H6l) is satisfied. 

In the log f vs. log n coordinates the formula (|44p yields the lines that resemble 
the behavior of the straight ones from the paper [8J. To exhibit it explicitly 
let us return to dimensionless variables and parameters and write 77 = 770 10 31 
where rj = 1 cP at which we typically have (A/D) ~ 10 (see previous Sec). 
Then we obtain 



log T = log 



u a u b \1_ /D\ 



2tt V 2vr \AJo 



log e 3x ( A\ log e 



where for the cubic (metastable) potential u a = 1, Ub = I and n = 6 while for 
the quartic (bistable) potential to a = a/2, u>b — 1 an d n = 4. The plots for the 
latter case obtained with the help of the formula (1471) are depicted in Fig. 5. 
However we should recognize that the lines obtained are in bad accordance 
with the results of the paper [8]. At realistic values A/D w 10 we obtain too 
large decrease of the reaction rate constant ~ 6 orders of magnitude instead of 
w 1 order of magnitude in [8j at the increase of solvent viscosity by 2 orders of 
magnitude. Besides the lines for this case deviate appreciably from the straight 
ones obtained in [8] . That is why in Appendix B we return to this problem and 
show that taking into account the corrections O(A) to the formula ( HO]) for 
the reaction rate enhancement improves the situation significantly and yields 
the results that are in good agreement with experimental data from the paper 



Now we consider a realistic cosolvent with finite molecular weight M. In this 
case both friction coefficients v and 7 are some unknown functions of solvent 
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viscosity. The cosolvent molecules with larger molecular weight worser pene- 
trate both into the enzyme active site where the reaction coordinate is located 
and into the enzyme body where the functionally important cranckshaft-like 
motion of favourably located peptide groups takes place. That is why with 
the increase of M these functions should grow more slowly the higher the M 
is. Thus cosolvent molecules with larger molecular weight should create less 
pronounced dependence on solvent viscosity than those with lesser molecular 
weight. In our opinion there are not enough experimental data to define these 
functions reliably on physical grounds. That is why we further resort to specu- 
lations of purely phenomenological character. We assume that these functions 
are of the form v oc r] S ^ M ' and 7 oc r/ a ( M ) with the requiremens that the in- 
decies of a power S(M) and a(M) obey \im.M^oS(M) = 1, liniAf^ocK^O = 1 
and are decaying functions with the increase of M. These requirements take 
into account the above mentioned features. 



Then the formula ff44l) yields 

f oc „,.., (4£ 



where 



(3{M) = 5{M) - '^P- (49) 



Regretfully the experimental data available for f3(M) do not allow us to distin- 
guish 5(M) and a(M) separately. We can only say that from the requirement 
/3(M) > we should have 5(M) > a(M)/A. The dependence of the parameter 
(3 on the values of 5 and a is presented in Fig. 6. One can distinguish two 
regions (see also Fig. 7): in the region / we have j3 > 0.75 while in the region 
77 we have (3 < 0.75. Accordingly three possible types of behavior are con- 
cievable. We call them normal, rare and exceptional and schematically depict 
in Fig. 7. Let us start from the case of large cosolvent molecular weight when 
both 6(M) and a(M) are small. For the rare type of behavior the decrease 
of M leads to a path in the a — 8 parameter space that is always left in the 
region /, i.e., the value f3(M) is always less than the limit value /3 max = 0.75. 
For the normal type of behavior the path crosses the fat line and enters the 
region II where (3(M) > 0.75. At the minimal realistic cosolvent molecular 
weight M = 62 for Ethylenglycol (see [8]) we have some a™ m and 8™ m at 
which we generally can obtain arbitrary 0.75 < /3 max < 1. Finally for the ex- 
ceptional type of behavior we can obtain a™ m so small and 5™ m so large that 
Pmax = 5™ in — a™ m /4 w 1, i.e., the approximate dependence F oc 1/r] emerges. 

Also the formula (1441 yields the dependence of the reaction rate constant on 
dielectric constant e. Such dependence has been known from the experiment 



for a long time pQ and is of the type In k cat = ... + const /e. The formula 
predicts In T — ... + const/e + (l/2)ln e. It is obvious that the latter term 
is practically unnoticeable on the background of the former one. That is why 
one can conclude that the formula (1441) yields correct dependence on dielectric 
constant in accordance with experimental data from [1] . 



8 Conclusions 



We suggest a dynamical mechanism that mediates the influence of solvent vis- 
cosity on the reaction coordinate at enzyme action. The mechanism is based 
on the fluctuations of the electric field in the enzyme active site. These fluctu- 
ations are produced by thermally equilibrium dynamics of protein structure, 
namely by rocking of the rigid plane of a favourably located and oriented 
peptide group (or may be a few of them). Such rocking causes the electric 
field of the dipole moment of the peptide group lying in its plane to undergo 
fluctuations in time. As the latter are thermally equilibrium they exist on 
the whole time scale of enzyme turnover and there is no problem to match 
them with the catalytic act. Namely the impossibility of survival for artificially 
constructed coherent oscillatory excitations in protein dynamics on the time 
scale of the catalytic act plagues many of the suggestions for dynamical con- 
tribution into enzyme action and causes justifiedly criticism of Warshel and 
coauthors [21], [32], [33]- Such excitations may really be created by energy 
released at substrate binding by an enzyme. They can exist for some life time 
in the form of, e.g., the so called discrete breathers either in protein highly 
regular secondary structure [35], [71] or in its whole irregular tertiary struc- 
ture [72]. However a protein is condensed media and any motion both in its 
interior and on its surface of thermally nonequilibrium characher must fade 
away rapidly (on the enzyme turnover time scale) due to dissipation. Long life 
times for discrete breathers obtained in [72] are due to unphysical assumption 
of friction for only surface elements of the three-dimensional network of oscil- 
lators. In our opinion the only possibility is to obtain some kind of influence 
on the reaction coordinate (having no counterpart for reactions in solution) 
from thermally equilibrium fluctuations. The structure of the protein enables 
peptide groups to implement high frequency thermally equilibrium rocking of 
their rigid planes (cranckshaft motion). This is a rotational motion and linear 
displacements of the atoms at it are negligibly small compared with both the 
size of the solvent molecules and the interatomic distances in protein. That 
is why this type of motion proceeds in the underdamped regime even in the 
condensed media of protein interior and mimics oscillations very much. 

The specific Fourier mode of the electric field fluctuations at some own fre- 
quency of rocking of the rigid plane of the peptide group (presumably u ~ 
100 cm' 1 ) resembles coherent oscillation and can be approximated by a har- 
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monic vibration. The latter is the RPV in our approach. The rocking of the 
rigid plane of the peptide group feels the microviscosity of the environment in 
the vicinity of the latter. This microviscosity is a function of solvent viscosity. 
Details of this function depend on the molecular weight of the cosolvent very 
much. For the hypothetical "ideal" cosolvent with infinitely small molecular 
weight (cosolvent molecules freely penetrate into enzyme and are distributed 
there homogeneously) the microviscosity is identical to solvent viscosity. For 
realistic cosolvent their relationship is more complicated and poorly known. 

The mechanism suggested in the present paper yields the required functional 
dependence ([1]) of an enzymatic reaction rate constant on solvent viscosity and 
the limit value (3 max ~ 0.75 that is in good agreement with the extrapolated 
one P m ax ~ 0.79 from the paper [8] (see Introduction). We have considered 
predominantly the limiting case of the hypothetical "ideal" cosolvent. To take 
into account the realistic molecular weight of the cosolvent one needs to know 
how its molecules are located in the enzyme both relative to the reaction 
coordinate and the rocking peptide group (or a few of them) interacting with 
the latter via electric field fluctuations. In this case both the friction coefficient 
for the reaction coordinate v and that for the rocking peptide group 7 become 
more complicated functions of solvent viscosity rj than simple proportionality. 
All the arguments from the papers pQ, [2], [15], [22], [23], [8], [17], [10] remain 
pertinent for this case. However we do not feel that there are enough data at 
present to tackle realistic cosolvent within our approach at a deeper level than 
purely phenomenological speculations. We are merely convinced that it should 
be done not within the Kramers' formula yielding k cat oc I/77 dependence but 
within a formula of the type of (jHJ) where the upper limit (3 max is already built- 
in, i.e., the dependence k cat oc \jrf with (3 < (3 max is obtained automatically. 

Regretfully at present we have to work with very limited tools from the 
Kramers' theory. The most natural formulation of our problem leads to a 
two noises model one of which is the oscillating noise (1231) . There are no re- 
sults for such model in the literature yet. Further development of the model 
will inevitably require its elaboration within the approach of the papers [62], 
[63] . At present stage we have to resort to artificial approximating the oscillat- 
ing noise by coherent oscillations. However even for the latter simplified case 
the theory of the Fokker-Planck equation (1381) does not provide us at present 
with a workable formula for the general set of the parameters A, D and Q. 
This fact does not enable us to evaluate reliably "how much can the electric 
field fluctuations in the enzyme active site contribute into the reaction rate 
enhancement ?" because small corrections in the exponent can lead to drastic 
overestimations. However we argue that the low frequency regime is relevant 
for enzymological problems. For this case the corrections 0(A) to the formula 
( HOI can be obtained explicitly (see Appendix B). They allow us conclude that 
the reaction rate enhancement by the suggested mechanism up to ~ 3 -j- 4 or- 
ders of magnitude is feasible. We argue that available results enable us to 
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consider the preexponent quite safely. As the functional dependence (DQ) of an 
enzymatic reaction rate constant on solvent viscosity manifests itself namely 
in the preexponent the employed formula (T4"0"j) seems to be sufficient for the 
purposes of the present paper. The discrepansies with experimental data are 
eliminated by taking into account the corrections 0(A) to this formula as the 
results from Appendix B testify. 

We conclude that the present approach grasps the main feature of the func- 
tional dependence ([1]) for enzymatic reaction rate constant on solvent viscosity 
and satisfactorily yields the upper limit for the fractional index of a power in 
it. 



9 Appendix A 



Here we give a brief estimate of the moment of inertia of the peptide group 
for rotation around the axis a that goes through the center of masses of the 
peptide group parallel to the bonds C l ~ x — C l and N l — C l a (see Fig.l). The 
peptide group is a rigid plane structure so that all atoms O, C, N and H 
lie in a plane. That is why the moment of inertia is simply the sum of their 
masses multiplied by the square of their distances to the axis a. The masses 
in atomic units (1 a.u.m. = 1.7 • 10~ 24 g) are tuq = 16, Vfic — 12, m^ = 14 
and irtH = 1. That is why the center of masses is shifted a little to the atoms 
C and O relative the middle of the bond C — N . The lengths of the bonds 
are loc — 1-24 A, 1 C n = 1-32 A and Inh — 1. A. That is why the distances 
to the axis a are approximately rp ~ 1.5 A, re ~ 0.3 A and r^ ~ 0.7 A and 
Th ~ 2 A. Substituting these values into the formula I = J^" 171 ^ we obtain 
I rj 7.6 • 10~ 39 g • cm 2 . Thus our rough estimate corroborate the precise value 
I w 7.34 • 10- 39 # • cm 2 from [52]. 



10 Appendix B 



We consider the low frequency regime Q << 1 that is argued in Sec. 6 to be 
relevant for our problem. For the particular case of moderately strong mod- 
ulation A « A c (where the bifurcational point A c has the values A c = 0.25 
for the case of the cubic (metastable) potential U(q) = q 2 /2 — q 3 /3 and 
A c = 2/(3v3) ~ 0.4 for the case of the quartic (bistable) potential U(q) = 
—q 2 f2 + q A /4:) the reaction rate enhancement is obtained in [73]. The formula 
obtained there is rather cumbersom and is not presented here to save room. 
The results obtained with its help for the case of quartic (bistable) poten- 
tial are depicted in Fig. 8. They testify that at our values D — 5 ■ 10~ 3 and 
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A/D ~ 10 (see Sec. 6) the reaction rate enhancement is ~ 3 -=-4 orders of mag- 
nitude. Thus the conclusion that the suggested mechanism is strong enough 
to affect enzyme action appreciably is justified. For the intermediate regime 
of moderately weak to moderately strong modulation D « A « \[D (our 
estimate A/D « 10 is within this range for practically important for enzymol- 
ogy range D « 3 4- 5 • 10~ 3 fj42l) ) the formula is simplified significantly and 
yields an explicit expression for the corrections 0(A) to the formula (J4"0l) . We 
denote 



P 



1 + — 




(50) 



where for the cubic (metastable) potential (CP) uj a = JU"(q a ) — 1, u> b = 
JU"(q b ) = 1 and q a = 0, q b = 1 while for the quartic (bistable) potential 
(QP) uj a = ^U"{q a ) = V2,uj b = ^U"(q b 



■li,i -= 1 and q a = -1, q b = 0. 
Then we have a simple formula for the reaction rate enhancenment 

15 {A 
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2*hrA 



- q a ) /2 - Ap 



=expi 



D 



q b -q a - Ap 



(51) 



where p is the constant given by (l50l . For the case of CP we have p ~ 1.632 
while for the QP we have p ~ 2.132. For both of them we have q b — q a = 1- 

From ( IBTT) we obtain a modification of the formula ( 1471) 



log T = log 
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+ log 
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(52) 



The results obtained with the help of this formula for the case of quartic 
(bistable) potential (uj a = a/2, cu b = 1 and n = 4) are depicted in Fig. 9. 
They are good agreement with the results of the paper [8\ because they yield 
straight lines for the decrease of the rate constants by approximately one order 
of magnitude at the increase of solvent viscosity by two orders of magnitude. 
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Fig. 1. Schematic picture of the peptide group. 
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Fig. 2. A look on the peptide group from the axis of rotation a explaining the 
definition of the angle <p (defined in ([3])). 
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Fig. 3. Model of the peptide group by oblate ellipsoid. 
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Fig. 4. Favourable orientation of the peptide group dipole moment p(t) relative to 
the reaction coordinate q. 
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Fig. 5. Dependence of the reaction rate constant f on solvent viscosity r? in the log T 
vs. log (77/770) coordinates given by formula (J47]) for the case of quartic (bistable) 
potential. The values of (A/D) from the down line to the upper one respectively 
are: 1, 6, 11, 16, 21. 
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Fig. 6. Dependence of the parameter (3 = 5 — a/4 (see (05J)) on the values of 5 and 
a. 
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Fig. 7. Three types of possible behavior of the parameters 5(M) and a(M) on 
cosolvent molecular weight M: n-normal, r-rare and e-exceptional. The fat line 
devides the parameter space to the region / where (3 = 5 — a/4 > 0.75 and the 
region II where (3 = 5 — a/4 < 0.75. 
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Fig. 8. The dependence of the escape rate enhancement on the driving amplitude to 
noise intensity ratio for the case of quartic (bistable) potential U(q) = —q 2 /2 + q/A. 
The values of the noise intensity D from the down line to the upper one respectively 
are: 5 • 10~ 3 ; 4 • 10~ 3 ; 3 • 10~ 3 ; 2 • 10~ 3 ; 1 • 1(T 3 . 
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Fig. 9. Dependence of the reaction rate constant f on solvent viscosity r] in the log T 
vs. log (77/770) coordinates given by formula fj52[) for the case of quartic (bistable) 
potential. The values of (A/D) from the down line to the upper one respectively 
are: 1, 6, 11, 16, 21. 
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